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ABSTRACT 



We propose a model for the gamma-ray binary LS 5039 in which the X-ray 
emission is due to the inverse Compton (IC) process instead of the synchrotron 
radiation. Although the synchrotron model has been discussed in previous stud- 
ies, it requires a strong magnetic field which leads to a severe suppression of the 
TeV gamma-ray flux in conflict with H.E.S.S. observations. In this paper, we 
calculate the IC emission by low energy electrons (7 e < 10 3 ) in the Thomson 
regime. We find that IC emission of the low energy electrons can explain the 
X-ray flux and spectrum observed with Suzaku if the minimum Lorentz factor 
of injected electrons 7 min is around 10 3 . In addition, we show that the Suzaku 
light curve is well reproduced if 7 min varies in proportion to the Fermi flux when 
the distribution function of injected electrons at higher energies is fixed. We 
conclude that the emission from LS 5039 is well explained by the model with 
the IC emission from electrons whose injection properties are dependent on the 
orbital phase. Since the X-ray flux is primarily determined by the total number 
of cooling electrons, this conclusion is rather robust, although some mismatches 
between the model and observations at the GeV band remain in the present 
formulation. 

Subject headings: binaries: close - radiation mechanisms: non-thermal - stars: 
individual (LS5039) - X-rays: binaries 
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Introduction 



Five gamma-ray bi naries have been id entified so far. T hey are LS 5039 (lAharonian et al. 



2005al ). LSI +61° 303 (lAlbert et a 



0632+057 (Aharonian et al 



2007 



2006h. PSR B125 9-63 (lAharonian et al 



Bongiorno et al 



2005bJ), HESS 



201lh . and 1FGL J1018.6-5856 



( jFermi LAT Collaboration et al.ll2012l ). These systems consist of a compact object, a 
neutron star or a black hole, and an OB star (LS 5039 and 1FGL J1018. 6-5856 have an 
O star while the others have a Be star). The nature of the compact objects is still open 
except for PSR B1259-63 whose primary star is known to be a young pulsar. 

LS 5039 consists of a compact object and a massive star whose spectral type is 06.5 
V an d whose mass is M* = 22. 91^9 M Q , and its binary period is 3.90 6 days (lCasares et al. 



20051 ). In addition, the system has the eccentricity, e = 0.24 ± 0.08 (jSarty et al. 



2011) 



and its separation changes from ~ 2.6 R* at periastron to ~ 4.2 R* at apastron, whe re 



20051 ). 



R* = 9. 31^0 Rq represents the radius of the companion O star (lCasares et al. 
Spectra and light curves of LS 5039 have been reported in the energy band of TeV gamma 
ray, GeV gamma ray, a nd X-ray, with H.E.S.S. array of atm ospheric Cherenko v telescopes 



(Ah aronian et al 



(ITakahashi et al. 



2006), Fermi Gamma-ray Space Telescope (lAbdo et al 



20091 ). and Suzaku 



20091 ). respectively. These observations show that the fluxes clearly 
modulate with its binary period, which implies that the high energy emissions are generated 
within the system. In addition, we can see from the observational data that the X-ray flux 
anti-correlates with the GeV gamma-ray flux and correlates with the TeV flux. We note 
further that the photon index in the X-ray band is ~ 1.5, which means that spectra of the 



electrons emitting X-rays have a power-law index ~ 2.0. 



Many authors have discussed the X-ray emission mechanism in LS 5039 



it is due to the synchrotron radiation ( Dermer fc Bottcher 

" II I I ' 



2006 



Zhang et al 



2009; 



2006 



Dubus 



2006 



assuming that 



Paredes et al. 



Takahashi et all 120091 : lYamaguchi fc Takaharall2010l : ICerutti et al 
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20101). Thev mode l ed the system in a 



Paredes et al. 



2006 



Zhang et al. 



2009 



microquasar scenario (IDermer fc Bottcherl 12006 



Takahashi et al 



2009 |), which in cludes jets from the 



20061 ). which includes a 



accreting compact object, or in a colliding wind scenario (jDubusI 
relativistic pulsar wind colliding with a stellar wind. In each scenario, they obtained results 
more or less consistent with the observations for all energy bands. 

However, in explaining the X-ray emission as the synchrotron radiation, there is a 
crucial pro blem. The magnetic field of around 3 Gauss is required to reproduce the Suzaku 



spectrum ( ITakahashi et al. 



2009 



Yamaguchi &: Takahara 



2010 



Cerutti et al 



20101), if the 



magnetic field is uniform and if a relativistic motion of a radiating medium is neglected. 
Under this magnetic field the cooling rate of the synchrotron process dominates that of the 
inverse Compton (IC) process for electrons with > 1 TeV. This leads to a severe suppression 
of the flux in TeV range, so that H.E.S.S. spectra cannot be reproduced, as shown in 



Cerutti et al 



(120 1Q[ ). 



On 
limited. 



;he other hand, the work t rying to explain X-ray emission as IC process has been 



Bosch-Ramon et al 



( 120051 ) suggested that the X-ray emission was due to the IC 



(and/or synchrotron) process from the jets and that the orbital variation of the accretion 
rate was responsible for the X-ray modulation before H.E.S.S. and Suzaku observations were 
reported. After this study, the IC model for the X-ray emission has not been discussed. In 
the light of the new observational data, we revisit the IC model to reproduce the spectra 
and the orbital modulation of the X-ray emission. In this paper we aim to find the condition 
to reproduce the Suzaku data when we take the same parameters as in our previous paper, 



Yamaguchi fc Takahara! ( 120101 ) (hereafter YT10). 



We describe the model in Section [2j and show obtained results in Section [3j where 
we demonstrate that the IC model can well reproduce the X-ray observations if certain 
conditions are satisfied. We discuss remaining problems and complexities in Section HJ 
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Finally, we summarize the study in Section [5j 



2. IC emission model for X-ray emission 

2.1. Model overview 

We assume that high energy electrons are accelerated with an isotropic distribution at 
the position of the compact object. Here, we do not specify the acceleration mechanism, 
i.e., the microquasar model or the wind collision model. Since the electrons are in the 
stellar radiation field, they lose energy by IC process. We neglect the synchrotron cooling 
by assuming that the magnetic energy density is much lower than the radiation energy 
density. This cooling process leads to modification of the electron distribution as given 
in the next section. On the other hand, the electrons emit X-rays and gamma rays by 
scattering off the stellar radiation. Since the stellar radiation field is anisotropic, the high 
energy emission is also anisotropic resulting from anisotropic IC scattering. This anisotropy 
affects the modulation of the emission. 

Here, we focus on low energy electrons whose IC photons are not subject to 77 
absorption, which have been neglected in previous works. In YT10, we injected electrons 
with 3 x 10 3 < 7 < 10 8 and calculated IC scattering, the 77 absorption and subsequent 
cascade process. In addition, we calculated the IC cooling down to 7 P = 3 x 10 3 . Thus, 



Takahashi et al. 



(120091 ) also 



we focused on the emission from higher energy electrons, 
concentrated on high energy IC emission (> 1 GeV). In this paper we calculate the emission 
from electrons with lower energy, which emit X-rays and MeV gamma rays through the IC 
process. 
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2.2. 



Electron distribution 



The energy distribution of injected electrons ni n j(7) is assumed as 



^inj 

(7) oc 7 p for 7 min < 7 < 7 



(1) 



where 7 m i n and 7 max are the minimum and maximum Lorentz factors of the injected 
electrons, respectively. Because gamma rays up to 30 TeV have been detected from LS 
5039, we assume 7 max = 10 8 . However, now we are interested in X-ray emission. Moreover, 
we have calculated GeV to TeV emission due to IC process by the Monte Carlo method and 
successfully reproduced the observations in YT10. Therefore, in this paper we calculate the 
IC emission from the electrons with 7 < 71 = 3 x 10 3 , which corresponds to the minimum 
Lorentz factor in YT10. In addition, we adopt the same power-law index of the injected 
electrons as in YT10, i.e., p = 2.5, and we assume that p and the coefficient in Equation ([T]) 
are constant in the whole orbit. 

The electrons severely lose energy due to the strong radiation field by the star. Given 
the energy density of the stellar radiation, the cooling time by IC process in the Thomson 
regime is 



where a represents the binary separation and 'peri' means periastron and 7 is the Lorentz 
factor of electrons. Equation (j2J) implies that injected electrons completely cool near the 
periastron in T orb /10 ~ 3 x 10 4 s, where T orb is the orbital period, while at the apastron 
electrons cool down to 7 C ~ 4 in T orb /10, where 7 C means the Lorentz factor of cooled 
electrons. 

Therefore, the electron distribution is modified by IC cooling. As a result, if 7 m i n > 7 C 




(2) 
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the distribution in the steady state is given by the following form, 

{7~ 2 for 7 C < 7 < 7 min 
(3) 
7 -(p+i) f or 7min < 7 < 7l 

The normalization of this distribution is determined by equating the number of electrons 
with 7! to that obtained at 71 by the Monte Carlo calculation in YT10. Therefore the 
electron distribution of Equation fl3]) is continuously joined to that of YT10. We note that 
11(7) in YT10 is just a boundary condition in this paper, so that the calculation below does 
not affect so much the GeV spectra. At the same time some mismatches between the model 
and observed spectra at the GeV band remain. 



2.3. Parameter 7 m i n 

We assume that 7 min is a free parameter so that we fit the IC emission spectra from 
these electrons with the Suzaku data. The X-ray flux observed by Suzaku anti-correlates 
with the Fermi flux. If we assume that 7 m i n is constant with the orbital phase, we expect 
that the fluxes at the two energy bands correlate with each other. This is because the 
orbital modulation is wholly due to anisotropic IC. In this case, the light curve of Fermi 
is reproduced by the anisotropy of IC process (see YT10), but that of Suzaku cannot be 
reproduced. Therefore, we should assume that 7 m ; n varies with the orbital phase so that the 
X-ray light curve fits the observed one. We note that this assumption is completely ad hoc 
at this stage, but it is possible that 7 min varies with the orbital phase because the condition 
of the electron injection could change with the orbital separation. 

Here, we discuss the relation between the X-ray flux and 7 m ; n in order to understand 
how we should set the variation of 7 m i n . We define F X)0 and n (7 x ) as the IC X-ray flux 
and the electron number spectrum at j x , respectively, expected when 7 min is constant with 
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the orbital phase, where 7 X represents the Lorentz factor of electrons emitting X-ray, i.e., 
7 X ~ 10. We note that F X;0 shows the orbital modulation similar to the Fermi flux -FceV 
because the modulation of F x correlates with that of the IC flux in the GeV band and 
because this GeV flux well reproduces the Fermi flux (see YT10). When 7 m ; n varies with 
the orbital phase, the expected X-ray flux F x and the electron number spectrum n(j x ) 
satisfy the relation 

^ = F Xi0 (4) 
^o(7x) 

Using Equation (J3J), we obtain 

n(^)=n(7l)f^ (^-) 2 ■ (5) 

V 7l / V7min/ 

Since we assume that the electron number at 71 varies in the same way as that in YT10, 
we set 

n(7i) = "o(7i)- ( 6 ) 

When 7 min is constant during the whole orbit, no (71) shows the same modulation as rio(7x), 
so that we obtain 

n(7i) oc no(7 x ). (7) 

Therefore, Equation (J5J) leads to 



1 mm 



n(7 x ) oc n (7 x )7 n 

As a result, Equation (jlj) reduces to 

F x oc F Xi0 7 mm 5 , (9) 

where we set p = 2.5. Therefore, we should increase 7 m i n when F x ,o(-Pbev) increases and 
vice versa, in order to reproduce the observed orbital modulation at X-rays. 
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2.4. Assumptions for the calculation 



We assume that the electrons which follow the energy distribution given by Equation 
([3]) are in a steady state for a given orbital phase because 7 X > 7 C . The stellar radiation 
is assumed to take a black-body distribution. The typical energy is 10 eV, so that IC 
scattering occurs in the Thomson regime for electrons 7 < 71. Therefore, we calculate 



nu merically the 



111 



Dubus et al. 



C sp ectra in the Thomson regime, which is done in the same manner as 



(120081 ). Here, we neglect the effect of secondary electrons created by the 



cascade process. 



Calculated Results 



In Figure [U the resultant spectra from low energy electrons (solid lines), the IC spectra 
from high energy electrons obtained in YT10 (dashed lines) and their sums (thin solid 
lines) are shown, where the spectra are averaged in the two orbital phases, around the 
inferior conjunction <p =0.45-0.9 (red lines and circles) and around the superior conjunction 
=0.9-0.45 (blue lines and circles). Here, <fi means the orbital phase, and = represents 
the periastron. The orbital variation of 7 m i n is given so that the light curve of the 
observation with Suzaku is reproduced (Figure Wj- We set the minimum Lorentz factor of 
cooling electrons as 7 C = 5, so that the cooling time of electrons with this Lorentz factor 
is shorter than a tenth of the orbital period at any orbital phase (see Section 12.21) . The 



inclination angle is taken a s % = 30° so t hat t 
parameters are taken from 



Sarty et al. 



lis model is consistent with YT10. The orbital 



(120111 ). The results do not change so m uch even if 



we take other allowable parameter sets (e.g. those given in 



Aragona et al. 



20091)- 



We determine the orbital variation of 7 m i n by fitting the calculated light curve to the 
observed one as shown in Figure [21 The comparison of the variation of 7 m i n with the Fermi 
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light curve is shown in Figure |3l where 7 m i n is shown with a solid line, _F Xj o is shown with a 
dashed line, and the flux obtained with Fermi is shown by crosses with dashed error bars. 
Figure [3] indicates that 7 min obtained from the fitting to the X-ray band shows almost the 
same variation as the IC flux with the constant injection, that is, 

7 min oc F Xi0 ~ F G ev- (10) 

In other words, if we assume the variability of 7 min as Equation (jTUl) . we can reproduce the 
observation with Suzaku. The reason why the fitting to the data of Suzaku yields such a 
simple relation is that the observed fluxes of X-ray and GeV bands are related as 

oc F"V (11) 

In addition to the orbital variation of 7 m i n , we note that this result gives a lower limit of 
7 min . If 7 min is lower than ~ 10 3 , calculated flux overpredicts the Suzaku one. 



4. Discussion 

We have shown in the previous section that the orbital modulation of the X-ray 
emission from LS 5039 is reproduced by the IC process if we give the ad hoc relation of 
Equation (fTUj) . However, it is not clear if such a dependence is expected. It is hard to 
assume that the minimum Lorentz factor of injected electrons varies in tandem with the line 
of sight. Nevertheless, since the existence of cooled electrons is inevitable, it is natural that 
the low energy electrons with 7 ~ 10 is responsible for the X-ray emission from LS 5039. 

In the model of this paper, the X-ray modulation is mainly due to the variation of 
the number of cooled electrons. Since the number of cooled electrons is determined by 
the injection rate J n iQ j('j)d'-f, the X-ray modulation is controlled by the injection rate. 
In this paper, we assume that p and the coefficient in Equation (CQ) are constant in the 
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whole orbit while 7 m ; n modulates. Here we note that the injection rate changes with the 
orbital motion. That is, the X-ray modulation is adjusted by making 7 min change in this 
paper. Alternatively, we may adjust the X-ray modulation by changing p and the coefficient 
because these two parameters are also related to the injection rate. 

While the Suzaku spectra are well fitted in Figure [TJ the calculated spectra in the GeV 
band match worse the observation since we adopt the YT10 data as higher energy spectra. 
The Fermi spectrum around the superior conjunction is now somewhat overpredicted 
even for the lowest two bins of Fermi spectrum. This means that the amplitude of the 
modulation in the GeV band is slightly overpredicted when the contribution from low 
energy electrons is considered. Therefore, in order to reproduce the Fermi observation we 
should set the smaller amplitude of the calculated light curve than the current model, which 
corresponds to a smaller inclination angle. Smaller inclination angle means a larger mass 
of the compact object, strengthening the suggestion in YT10 that it is a black hole. This 
modification of the inclination angle does not affect so much Equation f fTUj) although the 
amplitude of _F Xj o and thereby that of 7 m i n are slightly modified. We note that in this model 
TeV gamma-ray can be emitted by the IC process since this model is compatible with the 
model in YT10. 

Here, we mention the effect of the escape of electrons on the model. We assume in 
this paper that electrons cool down to 7 C without escaping. However, when taking into 
account the hydrodynamics of the system, e.g., the outflow of the jet in the microquasar 
scenario, if we assume that electrons are advected by the flow, they cannot cool radiatively 
so efficiently. For example, if electrons are advected at the light speed, the crossing time 
in the system t cross is ~ 100 s. Since the electrons cool during this crossing time, 7 C is 
determined by 

*ic,t = t cross j 

(12) 
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where £ic,t is given by Equation ([2]). Thus, we obtain 7 C ~ 300, which means that electrons 
do not emit X-rays through the IC scattering if advection is very fast. This issue may 
be solved in two ways. One is the adiabatic cooling where the spatial extent of electron 
confinement increases more than two orders of magnitude. The other is that electron 
injection itself continues to as low as 7 m i n ~ 10 with an index 2.0 below a break, then the 
photon index in the X-ray band may be reproduced. Since 7 C > 7 m i n , electrons emitting 
X-rays are not affected by the IC cooling. In these scenarios, the X-ray flux may modulate 
due to not only the anisotropic IC process but also the orbital variation of 7 C , but further 
study is required to reproduce the Suzaku observations. 

We comment shortly models based on the synchrotron emission. If we ascribe the 
X-ray emissi on to the syn c hrotro n radiation, the TeV observation cannot be reproduced 



as shown in 



Cerutti et al. 



( 120101 ). This discrepancy may be solved by a model in which 
the magnetic field is non-uniform. For example, X-rays are emitted in a compact region 
where the magnetic field is more than 3 Gauss and the electrons are populated in energies 
lower than 1 TeV, while TeV IC emission by the higher energy electrons takes place in 
an extended region with magnetic field strength as small as 0.1 Gauss. As an alternative 
scenario addressing this problem, we could consider a situation where the radiative matter 
moves along near the line of sight at relativistic velocity. In this scenario, the Suzaku 
spectra would be reproduced under the magnetic field less than 1 Gauss because of the 
relativistic beaming, and a severe suppression in the TeV band could be avoidable. 



Takahashi et al. 



( 120091 ) proposed the model in which the X-ray emission is explained 
by the synchrotron radiation. They assume that the adiabatic cooling is more efficient 
than the radiative cooling, so that the X-ray flux is proportional to the adiabatic cooling 
time. Thus, the X-ray light curve can be reproduced if one give properly the variation 
of the adiabatic cooling time. Moreover, since the adiabatic cooling does not change the 
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power-law index of the electron distribution, the spectral index in the X-ray band is also 
reproduced if one give 71^(7) oc 7~ 2 . In addition, they show that calculated TeV spectra 
match the H.E.S.S. observation. The suppression in the TeV band is avoided in this model 
because high energy electrons cool not radiatively but adiabatically, so that they are not 
affected by the synchrotron cooling. 

Many authors have constructed models in which the par ameters of inje c ted e lectrons 



vary with the orbi tal motion to reproduce the observations. 



Zhang et al. 



Paredes et al 



(|2006h and 



(120091 ) calculated the emission from the jet in the mi cro quasar scen a rio an d 



Paredes et al. 



()2006|) 



showed broad band spectra from radio to TeV gamma-ray band, 
assumed that the injection rate of the electron distribution varied with the orbital motion due 
to the change in the accretion rate, and investigated the orbital modul ation in each e nergy 



Zhang et al. 



J2OO9J) 



band. However, their calculation was based on the data before Suzaku. 
assumed that the electron power-law index and the accretion rate, which is equivalent to the 
injection rate, showed two values for two orbital phase, superior and inferior conjunctions, 
respectively. As a result, the Suzaku spectra are well reproduced but they d i d not m ention 



the orbital modulation for any 



energy band. 



Sierpowska-Bartosik fc Torres! (j2008al ) and 



Sierpowska-Bartosik fc Torres! (l2008bl ) calculated the IC scattering from electrons in the 
pulsar wind in the wind collision scenario and showed that their spectra and light curves 
reproduced the H.E.S.S. data. They assumed that the distribution of the injected electrons 
showed two different power-law indices for the two phase intervals, respectively, in order to 
reproduce the H.E.S.S. spectra. However, they have not investigated X-ray emission. 



Bednarek! (120111 ) assumed that in the wind collision scenario electrons were accelerated 
in two shocks separated by a contact discontinuity. He have shown that the IC emission 
from the pulsar-side shock and the star-side shock accounts for TeV emission and GeV 



emission, respectively. In his model the maximum energy of electrons varied with the orbital 
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motion, and the injection rate was constant. As a result, the normalization of the electron 
distribution is changed, so that X-ray emission by the synchrotron radiation modulates, 
although the flux varies inversely to the Suzaku data. 

The crucial difference between the IC model and the synchrotron model is in the MeV 
spectra. The IC model in this paper predicts that there is no dip in Me V spectra. In 



contra st, the s ynchrotron model pr edicts that there is a dip in MeV spectra (ITakahashi et al. 



()2009|), YT10, 



Cerutti et al 



( 120101 )). Therefore, future observations in the MeV band will 



judge which model is realized. 



5. Summary 

We assume that X-ray emission from LS 5039 is due to the IC scattering of electrons 
with 7 ~ 10. We calculate the anisotropic IC emission from electrons with the Lorentz factor 
5 < 7 < 3 x 10 3 in the Thomson regime. As a result, we have found that we can reproduce 
Suzaku flux if the minimum Lorentz factor of injected electrons 7 m ; n ~ 10 3 . Moreover, the 
X-ray light curve is well reproduced if a simple relation 7 min oc -F Ge v is satisfied, where -F Ge v 
represents the GeV flux. We have obtained these results based on YT10, where the index 
and the coefficient of the power-law distribution for injected electrons do not change with 
respect to the orbital phase. We may obtain different results when we adopt the model in 
which the index and the coefficient vary, but Suzaku data can be reproduced by adjusting 
these parameters, because the X-ray flux variation is primarily determined by the number 
of low energy electrons. Thus, we suggest that the X-ray emission obtained with Suzaku 
can be explained by the IC emission. However, we need more investigation for explaining 
the Fermi data. This model will be justified or rejected by future observations in the MeV 
range. 
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Fig. 1. — IC spectra by low energy electrons (7 < 3 x 10 3 ). The spectra are averaged in the 
two orbital phases, around inferior conjunction 0.45-0.9 (red lines and circles) and around 
superior conjunction 0.9-0.45 (blue li nes and circles). The observed spe ctra are shown b y 



circles with error bar for Suzaku (XIS, 



Takahashi et al. 



[20091) and Fermi (lAbdo et alJl2009h . 
IC spectra obtained by Monte Carlo method for 7 > 3 x 10 3 are also shown with dashed 
curves, which are equivalent to those in Figure 6 of YT10. We change 7 min as a function of 
the orbital phase (Figure so that the calculated X-ray light curve fits the observed one as 
shown in Figure [2j Thin solid lines show the sums of the spectra from low and high energy 
electrons in each phase. 
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Fig. 2. — Calculated (solid line) and observed (circles) X-ray light curves for the case as in 
Figure [TJ The vertical lines show the superior conjunction (supc) and inferior conjunction 
(infc). 
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Fig. 3. — Adopted phase dependence of 7 m i n (solid line) for the case of Figures [T] and |2j 
Dashed line shows -F x ,cb an d circles with dashed error bar show the observed GeV light 
curve. 



